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We demonstrate that quantum fluctuations can cause, under certain conditions, the dynamical 
instability of pure states that can result in their evolution into mixed states. It is shown that 
the degree and type of such an instability are controlled by the environment-induced anti-Hermitian 
parts of Hamiltonians. Using the quantum-statistical approach for non-Hermitian Hamiltonians and 
related non-linear master equation, we derive the equations that are necessary to study the stability 
properties of any model described by a non-Hermitian Hamiltonian. It turns out that the instability 
of pure states is not preassigned in the evolution equation but arises as the emergent phenomenon 
in its solutions. In order to illustrate the general formalism and different types of instability that 
may occur, we perform the local stability analysis of some exactly solvable two-state models, which 
can be used in the theories of open quantum-optical and spin systems. 

PACS numbers: 03.65.-w, 05.30.-d, 05.45.-a 


I. INTRODUCTION 


It is well-known that purity is exactly preserved dur¬ 
ing unitary evolution driven by Hermitian Hamiltonians. 
This property is natural for describing isolated quantum 
systems but in the case of open ones it is no longer com¬ 
pulsory Hi- From the viewpoint of the theory of open 
quantum systems, the isolated system is a mere theo¬ 
retical idealization - since in the real world all quantum 
wstems are embedded into a background of some kind 
[S, 1^ . Even the process of quantum measurement itself 
fits into this framework, because it involves interaction of 
the quantum system, which is being measured, with an 
external apparatus. Correspondingly, once the system is 
brought into interaction with its environment, such as a 
heat bath, dissipation usually increases its entropy and 
pure states are converted into mixed ones [1, @ . 

Recently, the dynamical behaviour of quantum purity 
and pure states has become of considerable research in¬ 
terest when studied within the framework of the non- 
Hermitian (NH) formalism 0,1^. Non-Hermitian Hamil¬ 
tonians find numerous applications in many areas of 
physics including studies of Feshbach resonances and 
decaying states, quantum transport and scattering by 
complex potentials, multiphoton ionization, free-electron 
lasers and optical resonators and waveguides. But the 
biggest area of application is the theory of open quan¬ 
tum systems where the anti-Hermitian part appears as 
a result of the interaction of systems with their environ¬ 
ment [ffil. In this paper we demonstrate the profound 
connection between NH Hamiltonians and dynamical sta¬ 
bility of pure states for the quantum systems of general 
type. We show that quantum fluctuations can cause the 
instability of pure states whose degree and type are gov¬ 
erned by the environment-induced anti-Hermitian terms 
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in Hamiltonians. 

The contents of this paper are as follows. In section 
m we provide a brief account of the density operator ap¬ 
proach for NH dynamics and formulate the essence of the 
stability problem for pure states. In section uni we de¬ 
rive the general equations which are needed for stability 
analysis of pure states for NH-driven systems and dis¬ 
cuss their generic features. In section HVl we study the 
linearized limit of the stability equations, which can be 
used for the analysis of local stability. In section |V1 we 
consider some two-level systems in order to illustrate the 
general formalism. Discussions and conclusions are given 
in section ED 


II. NON-HERMITIAN DYNAMICS 

If the Hamiltonian of a quantum system is a non- 
Hermitian operator, then it can be decomposed into its 
Hermitian and anti-Hermitian parts, respectively: 

H = H++H.=H+- it, (1) 

where H± = and f = f 1 is usually dubbed the de¬ 

cay operator. The probability-conserved time evolution 
of such a system is described by the normalized density 
operator p, which can be cast in the form 

/5 = f2/trH, (2) 

where Vi is called the non-normalized density operator. 
This operator is defined as a solution of the evolution 
equation 




tl 



( 3 ) 


where the square brackets denote the commutator and 
the curly ones denote the anti-commutator. This 
evolution equation can be directly derived from the 
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Schrodinger equation, see, for instance, ref. Q- Fur¬ 
ther, in this equation one can change from Cl to p, and 
obtain the evolution equation for the normalized density 
operator itself 

where we imply the standard definition for mean values, 
(A) = tr(/5i). 

From the mathematical point of view, equations ([3]) or 
together with the definition for computing mean val¬ 
ues, represent the map that allows to describe the time 
evolution of system © in terms of the matrix differen¬ 
tial equation which is defined on similar axiomatic foun¬ 
dations as the conventional master equations of the Li- 
ouville and Lindblad kind [H, [2^. According to this 
map, the Hermitian operator H+ = {H + takes 

over a role of the system’s Hamiltonian (cf. the com¬ 
mutator in equations (| 31 ) or ®) whereas the Hermitian 
operator F = i{H — H^)/2 induces the additional terms 
in the evolution equation that are supposed to account 
for new effects. In other words, a theory with the non- 
Hermitian Hamiltonian H is dual to a Liouvillian-type 
theory with the Hermitian Hamiltonian {H + H^)/2 but 
with the modified evolution equation, which thus be¬ 
comes the master equation of a special kind. This map¬ 
ping not only reveals new features of the dynamics driven 
by non-Hermitian Hamiltonians but also facilitates their 
application for open quantum systems (l^ . 

From the viewpoint of theory of open quantum sys¬ 
tems, the evolution equation for the non-normalized den¬ 
sity operator H effectively describes the original subsys¬ 
tem (with Hamiltonian i?+) and the effect of environ¬ 
ment (represented by F). Consequently, the evolution 
equation for the normalized density operator p effec¬ 
tively describes the original subsystem together with 
the effect of environment F and the probability flow be¬ 
tween the subsystem and reservoir. From the viewpoint 
of the subsystem alone, this flow is an essentially non- 
Hamiltonian process since it is described not by means 
of any kind of Hamiltonian but through the last term in 
the evolution equation Q . This makes NH models some¬ 
what similar to the Lindblad-type ones, where the effect 
of the environment is encoded in the evolution equation 
through the additional term often dubbed the dissipator. 
However, the important difference is that the Lindblad 
dissipator, which is a traceless operator by construction, 
does not affect the conservation of probability of the sys¬ 
tem whereas the last term in (|3]) restores the probability’s 
conservation which is otherwise broken by the anticom¬ 
mutator term. Besides, one can see that the last term 
in ® is nonlinear with respect to p, unlike its Lind¬ 
blad analogue. It is interesting that the appearance of 
nonlinearities in NH-related theories has been also sug¬ 
gested some time ago, although on different grounds of 
the Feshbach-Fano projection formalism [^ . 

Further, if one introduces the quantum purity V = 
tr(p^) then one can show that its time evolution is gov¬ 



erned by the equation 


|p = TO).l[{r)P-tr(Pf) 


( 5 ) 


where TZ{p) is the purity rate function Q. It is easy to 
see that the rate function vanishes identically in the case 
of Hermitian evolution (F = 0), but is otherwise an es¬ 
sentially non-trivial function of both the density operator 
and the anti-Hermitian part of Hamiltonian. As long as 
the density matrix pp of a pure state, which is defined via 
the idempotence p^ = pp, is the equilibrium point in both 
the Hermitian and non-Hermitian cases (he., TZ{pp) = 0), 
one might expect that any pure state is always preserved 
during the NH evolution. 

However, this could only be possible if one disregards 
the other important player in the quantum realm - quan¬ 
tum fluctuations. Indeed, if an (initially) pure state 
is not protected against the fluctuations that can alter 
its purity, then during time evolution it will be driven 
away from being pure, no matter how small these fluc¬ 
tuations initially were. This phenomenon of dynami¬ 
cal instability is not directly seen in the evolution equa¬ 
tion for the density operator, but emerges via solutions 
thereof. It is somewhat analogous to the spontaneous 
symmetry breaking in field theory - except that here one 
deals not with the actual field potential, but with some 
kind of the fictitious-particle potential function (usually 
dubbed as the Lyapunov function candidate) that deter¬ 
mines whether the equilibrium point pp is stable or not. 
This also makes NH models different from the Lindblad- 
type ones, in which the non-conservation of purity and 
pure states directly follows from the underlying evolu¬ 
tion equation. Besides, the Lindblad master equation 
approach has a different range of applicability because it 
implies a number of certain approximations, such as the 
Markovian, Born and rotating-wave ones, whereas NH 
Hamiltonians often appear in theories in a more direct 
way, one example to be the Feshbach-Fano projections 

m 

It is worth mentioning also that the above-mentioned 
kind of chaotic behaviour should not be confused with 
the notion of “quantum chaos”, which is currently being 
reserved in the field of research devoted to how classical 
chaotic dynamical systems can be described by means of 
methods and concepts of quantum mechanics (26l - [^ . 


III. GENERAL STABILITY APPROACH 

In order to develop the general approach for perform¬ 
ing the stability analysis of pure states, let us introduce 
the non-purity (mixedness) operator 

M = p-4, (6) 

whose trace is known as the linear entropy 
Sl = trM =l-V. 


( 7 ) 
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In order to consider variations of the density operator 
around some pure state /5p, we perform the decomposition 

P = /5p + A, (8) 

where A is the variation operator. As long as the main 
properties of the density matrix should be left intact by 
such decomposition, A must be Hermitian and traceless; 
this automatically ensures the well-defined probability 


and real mean values of operators. One obtains that 

p^ = Pv + {pp, A} + A2, M = A - {pp, A} - A2, (9) 

such that the operator M is clearly a measure of deviation 
of a state from being pure. Consequently, equations (|3]) 
and yield the equations that are more suitable for 
stability analysis. 


dt 

dt 

d 


-A 


M 


dt 


Sl 


i 

h 

H+,A 

- |f, A| -1 - -Pp(F)a + ^ ((F)p -1 - (F)a) a. 

i 

h 

'h+,m 


(Pp + A)A(^) + A^r) 


((r)p + (r)A)5L-tr(fM) 


( 10 ) 

( 11 ) 

( 12 ) 


where (r)p = tr(ppf) is the average of the operator F 
with respect to the unperturbed density operator pp, 
(Ha = tr(Ar) is the difference between the averages 
(r) and (r)p, and we have used the notation = 
FA — tr(rA)J with I being the identity operator. One 
can easily verify that both the tracelessness and hermitic- 
ity of the variation operator A are preserved during time 
evolution. 

Equations (nnii-diii) are matrix differential equations 
that can all be used for the stability study, but it must 
be emphasized that equations m and uni) describe only 
those variations that can potentially lead to the transi¬ 
tion of a pure state pp into a mixed one (dubbed as the 
“mixing” fluctuations in what follows), whereas equation 
alone governs the variations of a general type, re¬ 
gardless on whether they alter the purity of pp or not. 
Therefore, the equations m and (HU) will be of spe¬ 
cial interest here. It is easy to see that one can apply to 
them the standard methods of dynamical stability analy¬ 
sis such as the Lyapunov or Vakhitov-Kolokolov criteria 

dg. 

However, some generic features can be noticed straight 
away. In particular, equation dm reveals that stabil¬ 
ity against the “mixing” fluctuations essentially depends 
on the result of competition between the terms (F) and 
tr(fM)/S'L = ((F) — tr(p^f))/S' l. If their difference is 
non-negative at any time, then during evolution the fluc¬ 
tuations will keep the initial state away from being pure. 
If the difference is negative then the fluctuations will get 
suppressed. 

Yet another property, which can be immediately seen 
from equation m, is that, in the absence of the anti- 
Hermitian part of the Hamiltonian, random fluctuations 
of a given pure state would never get suppressed but their 
magnitude remains at much the same level. This means 
that in order to fully suppress fluctuations, suitably cho¬ 


sen anti-Hermitian terms must be present to the Hamil¬ 
tonian. In those cases the anti-Hermitian part would en¬ 
sure the full stability of pure states even if it is negligibly 
small compared to the Hermitian part. 


IV. LOCAL STABILITY 

It is clear that the local instability of a pure state 
against fluctuations leading to the mixing of a state does 
yet not imply global instability (i. e., when the purity 
of a perturbed pure state never goes back to its origi¬ 
nal value I). However, this case is still interesting from 
a physical point of view: quantum systems might exist 
in which the purification time {i. e., the time of return 
back into a pure state) can be larger than the lifetime of 
the system itself or the ultimately possible time of mea¬ 
surement/observation of the system. Yet another appli¬ 
cation area of local instability is that it can point to the 
presence of singular points at which the density matrix 
components diverge; thus indicating that the underlying 
system becomes critically unstable. In any of these cases, 
local instability of a pure state might become the domi¬ 
nating chaotic-type phenomenon, which can strongly af¬ 
fect the physical properties of systems’ evolution. From 
the technical point of view, the advantage of the study 
of local instability is that one does not have to know the 
whole solution for the density operator, but only a few 
characteristic exponents. 

The analysis of local stability for a given state of the 
NH quantum system can be made based on linearized 
equations, which are easier to solve or analyze using well- 
known methods from the theory of stability and dynam¬ 
ical chaos. Indeed, once we assume that the variation 
operator is small, A = Sp, we can perform the lineariza¬ 
tion with respect to it. Then equations (nnn-dni) yield. 
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respectively: 


J 


^ A- 

^SM 

dt 


i 

h L 

i 

h L 


~I (r)p(5p + pptr(^pf) +0{Sp'^), 

H+,SM^ + l [2{T)pSM + ppA^^'^ + p^pSp\ +OiSp^) 


h i 


{T)p SSl - tr(f5M) + OiSp"^) 


(13) 

(14) 

(15) 


where SM = = 6p — {pp,6p} and 6Sl = 

tr{5M) = —2ii{pp 6p) = —6V are the variations of the 
non-purity operator and its trace, respectively; O de¬ 
notes the terms that can be neglected in the linear order 
of approximation. 

As discussed above, in order to study the stability 
of a system against the “mixing” fluctuations, one has 
to consider primarily equations m cinci (|15p . In the 
leading approximation da can be rewritten in the form 
— AX, where X is a column vector constructed 
out of the linearly independent components of the ma¬ 
trix 6M, and A = A{pp,H^,r,t) is a characteristic 
matrix. If all real parts of its eigenvalues are negative 
then the pure state pp is locally stable against small 
“mixing” fluctuations. This condition is equivalent to 
the matrix A^P -|- PA being negative definite for some 
positive definite matrix P = P^. The corresponding 
Lyapunov function candidate can be determined then as 
V (X) = (X)^P X. Further classification of instability 
types (“center”, “node”, “saddle”, “spiral”, etc.) can be 
done for a specific system by analysis of the characteristic 
matrix’s eigenvalues [s^. 

V. EXAMPLE: TWO-LEVEL SYSTEMS 

In order to illustrate the formalism, let us consider 
some models which exhibit, under certain conditions, the 
dynamical instability of pure states. 

The NH two-level systems (TLS) are simple yet very 
useful physical models (3lj | , which can provide a clear vi¬ 
sualization of different kinds of stability that may occur. 
One can check that for two-level systems the non-purity 
matrix has only one independent component and can be 
written in the form 

M = ^Sj, (16) 

then one obtains 

SM = ^SSj, (17) 

i. e., the scalar SSl provides complete information about 
the local stability of a given pure state pp, whereas the 


characteristic matrix A can be reduced to a 1 x 1 matrix. 
The linearized equation m takes the simple form 

jL|^S. = Afe.f).|(2(r),-trf). (18) 

where A = A{pp, T) is the characteristic exponent. Ac¬ 
cording to the approach, the state pp is locally stable 
against the “mixing” fluctuations if A{pp, T) < 0, and 
locally unstable otherwise. To illustrate this, let us con¬ 
sider the following exactly-solvable models. 

A. Tunneling models with non-Hermitian detuning 

This is a set of NH models whose Hermitian and anti- 
Hermitian parts of Hamiltonian are, respectively: 

H+ = —hujdx, r = hX&z, (19) 

where d’s are Pauli matrices, positive parameter oj is re¬ 
lated to the matrix element for tunneling between two 
wells, and A is a constant parameter, real-valued but oth¬ 
erwise free, that can be viewed as the imaginary coun¬ 
terpart of the detuning parameter [s^. Such models are 
popular in many areas of quantum physics, including the 
theory of open quantum-optical systems (in particular, 
when describing the direct photodetection of a driven 
two-level atom interacting with the electromagnetic field 
- in which case A would be related to the atomic damp¬ 
ing rate) iEl and non-Hermitian quantum mechanics 
with real energy spectra [H, . 

Here we would like to determine for which values of 
the parameters lu and A a pure state, say, the one defined 
by 

pW = |e){e| = , (20) 

is locally stable against the “mixing” fluctuations. In 
case of quantum-optical or spin systems, this information 
could be instrumental for determining whether a system 
can undergo, respectively, the spontaneous emission or 
spin-flip transition. 

The characteristic exponent A, which can be immedi¬ 
ately computed from equation (UHl), turns out to be equal 
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to A, up to a positive factor. Therefore, we expect the 
state p'p'^ to be locally stable against small fluctuations 
for the models with negative A = A/w, and locally unsta¬ 
ble otherwise. 

To verify this, one needs to solve the evolution equation 
assuming the perturbed initial conditions: 

p(0)=4i)+,5p(o), (21) 

where, according to the approach, the variation matrix 
can be chosen in the form 


= 4 ), ( 22 ) 


with 5iS being arbitrary real-valued numbers (for sim¬ 
plicity we have omitted the off-diagonal imaginary com¬ 
ponents of A). The positivity of the total density oper¬ 
ator p(0) can be always ensured by imposing constraints 
on the values of variation matrix’s components. However, 
here we assume the fluctuations to be arbitrary enough 
in a sense that we do not postulate them to exactly pre¬ 
serve the positivity of the density operator at the initial 
moment of time, similarly to the approach [l^ . Instead, 
one might be interested to see whether this property can 
be dynamically restored during evolution if it has been 
initially broken by fluctuations. 

Solving the evolution equation (jT]) with the initial con¬ 
dition (EH), we obtain the following expression for the 
normalized density operator: 


p{t) = WWTvCTx 


fyit) 


fzit) 


' F{t) 2F{t) 
where we have denoted: 

fx{t) = 


(23) 


(24) 


fy{t) = sinh (fir) ^p 2 cosh (/rr) — Asinh (/rr) , (25) 
fz{t)=P /ip 2 cosh (2^r) — A sinh (2^r) , (26) 

F{t) = A^cosh(2/rT) — p 2 AtAsinh (2/rT) — 1, (27) 


where t = ojt, X = X/uj, p 2 = 262 + 1 and fi = ■\/A^ — 1. 

The evolution of purity for this solution is presented 
in Fig. [TJ The top (bottom) row of the figure shows 
models for which the initial state is unstable (sta¬ 
ble) locally, as predicted by theory. Indeed, one can see 
that for positive values of A, which has the same sign as 
A, the magnitude of the linear entropy initially increases 
from its initial value, whereas for negative A’s the mag¬ 
nitude of the linear entropy decreases straight from the 
beginning. We reiterate here that the local stability or 
instability does yet not imply the global one. However, if 
the measurement time or the lifetime of a system is lim¬ 
ited and smaller than the purification time (which is of 
the scale uj~^ here), then an observer would not be able 
to empirically distinguish between the local and global 
(asymptotical) types of instability. 



FIG. 1: Linear entropy Sl versus time (in units of uj '^) for 
the family of models m, at different values of A and initial 
perturbations of the state p^K The value of 5i is 0.01 for 
all curves, the values of S 2 are: —0.02 (solid curves), —0.01 
(dashed curves), 0 (dash-dotted curves), 0.01 (dotted curves) 
and 0.02 (dash-double-dotted curves). 


Further, the left (right) column of Fig. [T] shows mod¬ 
els for which the initial state p^^ is unstable (stable) 
asymptotically. Asymptotical stability cannot be as¬ 
sessed based on the linearized equation (fTSl) but requires 
the usage of a general approach dinD-diii), which is quite 
bulky and thus would bring us beyond our purposes here. 
Instead, in Figs. [Hd and [TJi we demonstrate that cor¬ 
responding models contain the mechanism that asymp¬ 
totically suppresses quantum fluctuations, whereas Figs. 
[Hi and IH reveal that in some models fluctuations may 
keep the state away from being pure in an oscillatory 
way. As long as both the conventional von Neumann en¬ 
tropy and the NH-adapted entropy (defined in Ref. mi) 
would also oscillate in those cases, one can regard this be¬ 
haviour as the quantum non-Hermitian analogue of the 
self-organization phenomenon. 


B. Tunneling models with analytically continued 
matrix element 

This is a family of models where 

iJ+ = -hioax, F = hpa^, (28) 

with r] being a constant parameter, real-valued but oth¬ 
erwise free. Such models describe physical cases when a 
quantum environment of some kind (external fields or 
vacuum oscillations) effectively shifts the value of the 
tunneling parameter into a complex domain: ui ui+ip. 
Here we would like to determine for which values of 
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(a) »? = -2 



(b) 7) = 2 


FIG. 2: Linear entropy Sl versus time (in units of uj for 
the family of models (I28II . at different values of fi and initial 

( 2 ) 

perturbations of the state pp . The value of 82 is 0.01 for 
all curves, the values of are: —0.02 (solid curves), —0.01 
(dashed curves), 0 (dash-dotted curves), 0.01 (dotted curves) 
and 0.02 (dash-double-dotted curves). 


the parameters u) and 77 the pure state 




(29) 


is locally stable against the “mixing” fluctuations. From 
the physical point of view, it might be used to determine 
whether a given quantum superposition of states is pro¬ 
tected against the spontaneous decay into a mixed state. 

The characteristic exponent (1181) . turns out to be equal 
to 77 , up to a positive factor, hence, we expect state pp to 
be locally stable against small fluctuations in the models 
with negative fj = rj/uj, and locally unstable otherwise. 

To verify this, one needs to solve the evolution equation 
(HI assuming the perturbed initial conditions 


m = + (30) 


where 5p(o) is given by (0^ . We obtain the following 
expression for the normalized density operator: 


-/.N 9x{t) 

pit) = -zrTTT-TCr: 


9vit) ^ , 9z{t) 


2G(t) G{t) 


G{t) 


(31) 


where we have denoted: 


9 x{t) = Pi cosh ( 277 T) — sinh {2fjT), (32) 

gy{t) = 62 sin {2t), (33) 

gz{t) = 62 cos {2 t), (34) 

G{t) = cosh ( 2 ? 7 r) — pi sinh ( 2 ? 7 r), (35) 


where r = wt, 77 = 77 / 0 ;, and pi = 2Si + 1 . 

The evolution of purity for this solution is presented 

in Fig. HI Figured^ illustrates the model for which the 
( 2 ) 

initial state pp is stable locally, as predicted by the lin¬ 
earized approach. It also turns out to be stable globally, 
which makes this case similar to the one discussed earlier, 
regarding Fig. [TJi. 

In Fig. Ht* we show the model for which some flue- 
tuation modes of the initial state pp grow indefinitely, 
so that the system evolves into a singularity. This is 


another kind of dynamical instability, which is different 
from those discussed in the previous set of models. It 
illustrates that for some models, quantum fluctuations 
either get significantly amplified during certain interval 
of time (solid and dashed curves in Fig. H};) or destabi¬ 
lize the original state up to the full destruction (dotted, 
dash-dotted and dash-double-dotted curves). 


VI. CONCLUSION 


Using the density operator approach for non-Hermitian 
Hamiltonians, we have demonstrated that quantum fluc¬ 
tuations can cause, under certain conditions, the instabil¬ 
ity of pure states, which is controlled by the environment- 
induced anti-Hermitian parts of Hamiltonians. This is 
drastically different from the Hermitian case where both 
the purity’s value and pure states are preserved during 
time evolution. 

It is shown that the instability of pure states is not 
preassigned in the evolution equation but arises as the 
emergent phenomenon in its solutions. We have derived 
the equations that are necessary to study the stability 
issues of any quantum system, regardless of the number 
of its degrees of freedom, dimensionality of Hilbert space, 
etc. Thus, the formalism and main results are applicable 
for systems described by non-Hermitian Hamiltonians of 
general type. 

Finally, in order to illustrate the different types of in¬ 
stability that may occur, we have considered some ex¬ 
actly solvable two-state models. Apart from being in¬ 
structive examples on their own, these models can be 
used in a theory of open quantum-optical and spin sys¬ 
tems where our stability analysis might be helpful in 
achieving a better understanding of such phenomena as 
the spontaneous emission, decay or spin flip (which, as a 
matter of fact, could be one of the directions for future 
work). By means of those models we have visualized the 
main types of stability of a pure state against small “mix¬ 
ing” fluctuations that may occur in the open quantum 
system: fluctuations get permanently suppressed with 
time (the state is locally and globally stable), fluctua¬ 
tions amplify during a finite period of time but eventually 
get suppressed (the state is locally unstable but asymp¬ 
totically stable), fluctuations never get suppressed with 
time but stay bound by a finite value (the state is glob¬ 
ally unstable), and, finally, fluctuations indefinitely grow 
thus leading to the singularity and critical instability of 
the system (the state and system are globally unstable). 
It should be also noticed that for the two-state models 
studied above the fluctuations, which cause the purity to 
acquire unphysical values, either get suppressed or lead 
to the overall instability of the system. 
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